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Abstract 

The convergence problem for scattering states is studied in detail within the framework of the Algebraic 
Model, a representation of the Schrodinger equation in an L 2 basis. The dynamical equations of this model 
are reformulated featuring new "Dynamical Coefficients", which explicitly reveal the potential effects. A 
general analysis of the Dynamical Coefficients leads to an optimal basis yielding well converging, precise and 
stable results. A set of strategies for solving the equations for non-optimal bases is formulated based on the 
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, asymptotic behaviour of the Dynamical Coefficients. These strategies are shown to provide a dramatically 

improved convergence of the solutions. 
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In the quest for solving the Schrodinger equation for both bound and continuum states, square-integrable bases 
have been repeatedly used. For bound states this turns out to be a traditional way to obtain the spectral 
properties of quantum systems. It has been shown however that a single representation of the Schrodinger 
| equation in terms of an L 2 basis can be formulated that allows for a description of both bound and continuum 
states §|,|,|,|. 

A version of such a formulation is called the Algebraic Model of the Resonating Group Method (hereafter 
referred to as the Algebraic Model or AM). It was originally tailored to treat clusterized problems, but can 
be used for all kinds of quantum mechanical many-particle configurations without major modifications. The 
AM has been formulated in terms of different types of bases, depending on the more specific features of the 
quantum system considered. One very important feature of the AM is the fact that the boundary conditions 
of the system are translated from a coordinate space context to the context of expansion coefficients, and are 
explicitly incorporated in the dynamical equations 0, |[ |[ ||. 

We will consider a specific AM formulation (i.e. a specific L 2 basis) to elucidate the analysis of the method. 
The methodology used for this analysis is of a general nature however, and can be repeated for other bases. 
The specific AM version chosen in this paper features an oscillator basis, and is in particular very suitable 
for obtaining both the bound and continuum spectra of nuclear systems with very different configurational 
properties; where appropriate, nuclear spectroscopic units will therefore be used. The choice of this AM version 
is mainly due to the background of the authors. Results of the AM approach considered here have already 



1 



been reported on |]6| [7J , and show in particular how the coupling of cluster and collective configurations can be 
treated seamlessly in such a description. 

Where necessary, a specific form for the potential operator will be used. We will consider a Gaussian 
form in this work. Although this again is a popular potential form in nuclear spectroscopic calculations, it 
is also an interesting functional form to approximate a large variety of potentials by discrete and continuous 
superpositions. 

We will concentrate our attention to the solutions of the dynamical equations for scattering situations only, 
as these are much more involved than the bound state problems. Indeed, the latter can be well approximated 
by a simple diagonalization of the energy matrix, as is well known. 

We will specifically discuss how strongly the convergence of the solutions of an AM system depend on the 
parameters of the problem. In particular the dependence on the precision and the convergence properties of the 
solutions of both the oscillator radius of the basis (the parameter unambiguously fixing the square-integrable 
basis) and the form of the potential energy contribution, will be treated in detail. 

The convergence is crucial for obtaining stable approximations to the solutions of problems expressed in 
terms of an (infinite) set of basis functions, in our case an L 2 basis. This problem was repeatedly investigated, 
mainly for bound-state solutions. As applications to scattering problems, expressed in an L 2 basis, appeared, 
several algorithms were suggested to accelerate the convergence of the results within a restricted subset of the 
basis. For instance, Heller and Yamani || used "Kato correction" , and Revai et al. J9| introduced the "Lanczos 
factor" . A more intuitive approximation was proposed by Fillipov et al. |(| . 

The analysis of the AM equations presented in this work will be shown to lead to (1) an optimal choice for 
the basis, given the potential, yielding well-converging and stable solutions of the AM system, or, if the optimal 
basis cannot be used, to (2) algorithmic procedures for solving the AM system in an acceptable and controlled 
approximation. These algorithmic procedures will depend on the specific choice made for the basis (i.e. the 
oscillator radius chosen), but also on information on the asymptotic behaviour of the expansion coefficients, 
as well as on the dynamical equations themselves. In this way the algorithm used will depend on the physical 
properties of the system considered. 

In a forthcoming paper the different strategies for solving AM equations will be applied explicitly to a 
number of problems from nuclear, atomic and molecular physics. 



2 The AM in an Oscillator Basis 

Choosing an oscillator basis to describe some specific Hilbert subspace in which to solve the Schrodinger equation 
leads to the following form of the latter in terms of a system of linear equations: 

oo 

<n \ H — E \ m> c rn = (1) 

m— 

where the coefficients c m are the expansion coefficients in the oscillator basis of the wave function corresponding 
to the energy E 

oo 
n=0 

and subject to a typical boundary condition. For simplicity we have omitted the angular momentum quantum 
number as well as the energy dependence of the coefficients c n . 

The boundary conditions of quantum systems are traditionally expressed in coordinate space, but can also 
be formulated in terms of the expansion coefficients c n . Indeed, for very large n, the dynamical equations are 
reduced to a simple (in an oscillator basis a three-term recurrence) form containing the kinetic energy operator 
solely. The equations can therefore be solved analytically for very large n [|[ f|, ||. The above mentioned 
reduction of the equations only involves the supposition that the potential energy matrix elements vanish for 
very large n, which has been shown to be an acceptable approximation for relatively short-range interactions. 
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One obtains the following asymptotic behaviour for bound states: 

<4°»> ~ v 7 ^ exp(~KR n )/R n (3) 
(k = \J-2mE/h 2 ) 

and for continuum states: 

<4°*> ~ v /^(i i (fcE n )+tan(5) n L {kRn)) (4) 

(Jfe = \j2mE/h 2 ) 

where j'l and ri£ are the traditional Bessel and Neumann special functions. 

A striking resemblance with the asymptotic forms of a wave function in coordinate representation is observed, 
by replacing the radial coordinate in the latter by the discrete value R n = ^/An + 2L + 3. A heuristic argument 
for this observation is that R n corresponds to the turning point of the oscillator in state | n > with angular 
momentum L. 



3 A Simple Scheme for Solving the AM System of Equations 

In the previous section we introduced the infinite dimensional system of linear equations obtained from the AM 
formulation, to be solved subject to a proper boundary condition. As indicated earlier we will concentrate on 
the scattering situation, and therefore only consider the asymptotic behaviour of the c„ corresponding to the 
continuum boundary condition. 

A scheme for solving the linear system of equations quite naturally presents itself. Under the assumption 
that the matrix elements potential energy < i | V \ j > vanish for sufficiently large values of one of the basis 
state indices i or j, one chooses a limiting value N to set this vanishing point in terms of the basis states. In 
the region where the potential matrix elements are neglected, the expansion coefficients c n are then given by 
(||), and can be written as 

c„ = c<+) + tan(<J) (n > N) (5) 

assuming N to be sufficiently large, with 

<4 +) = \fR^3L{kR n ) (n>N) 

(£-> = ^R~n L (kR n ) (n>N) (6) 
The choice for N divides the linear system in three different regions: 

• a finite number of N equations with < n < N, in which the potential energy matrix elements are fully 
taken into account. These equations correspond to the "internal region" in terms of the basis states 

• an infinite number of equations corresponding to n > N, in which the potential energy matrix elements 
are neglected. These equations correspond to the "asymptotic region" , and are trivially fulfilled due to 
the boundary condition 

• the equation with n — N, in which the potential matrix elements arc neglected. This equation corresponds 
to the "matching condition", as it couples the internal region through coefficient cjv-i with the asymptotic 
region through the phase-shift S 

This scheme amounts to solving the following N + 1 dimensional system of linear equations for the N 
coefficients c„ with n = 0, 1, .., N — 1, and the phase-shift 5: 

/ H 00 — E ■ ■ ■ H 0<N _i \ / c \ 

V o 
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where 



T, 



N 



N 



(Tn,n — E) + Tjv,jv+i Cx +1 
{T n ,n ~ E) + T NjN+1 c^|j 



(8) 



and T stands for the kinetic energy operator. 

As the potential matrix elements do not actually drop to zero exactly for n > N , one should vary the value 
of N to test the stability of the solution. It turns out that this stability strongly depends on the specific problem 
considered. As an example, the solution of the linear system of equations for a nuclear two-cluster problem, in 
which the distance of the clusters is considered as the degree of freedom, shows a rapid convergence in terms of 
N. A solution for a monopole description of the nucleus, in which the radius of the nucleus is considered to be 
the prominent degree of freedom, shows a very slow convergence in terms of N. These results indicate that one 
should be very careful when omitting potential energy matrix elements, and that a proper study of the form of 
the equations is necessary. 



4 An Analysis of the AM Equations 
4.1 A Reformulation of the AM Equations 

To study the properties of the (in principle infinite dimensional) linear system (Q) to be solved, we will rewrite 
this set of equations using the following substitution for the expansion coefficients c n : 

Cn = C n + C n (9) 

By this substitution, the coefficients c„ are considered to represent a deviation from the asymptotic behaviour, 
i.e. the coefficients ci as ' . The coefficient Cn^ then quantifies this deviation, which in particular will be zero in the 
true asymptotic region (i.e. for very large n). The first term in (^) is responsible for the long-range behaviour 
of the system. The second term corresponds to the short-range correction caused by the potential; in other 
words, in coordinate representation the coefficients {ci°^} would represent that part of the wave- function that 
is dominated by, and within the range of, the potential. The knowledge of oh^ as a function of n thus provides 
a key element for determining a proper index N distinguishing the internal from the asymptotic region in terms 
of the basis functions. 

Rewriting the original AM linear system of equations (Q) in the unknowns ({c n }), yields an equivalent 
linear system in the unknowns ({cn°^}, 5), where 5 is defined by the following form of the asymptotic expansion 
coefficients c„ 

c^)= C W+tan(5)c(-) (10) 

Because the asymptotic coefficients now appear for all n in this representation, they should be properly 
defined. In order to do so, we consider the coordinate representation of the outgoing asymptotic wave- functions 
and which are originally defined as the two linearly independent solutions of 

(f-J5)* = (11) 

where T is the kinetic energy operator. "I'M is commonly called the "regular" solution, and behaves properly 
for all r. is the "irregular" solution, and has an irregular (infinite) behaviour near the origin r = 0. To 
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provide a regular character at the origin for both and \ we redefine their equations in the following 
way, as was suggested earlier by Heller and Yamani 

(T = 

(f = /3 $ (12) 

The non-zero right-hand side in Jl^), in which $o( r ) =< r | > represents the zero-quanta oscillator state, and 
where /3q equals 

exp(fc 2 /2) / 2 1 

ft-*" fc L + i y r (L + i/2) r(-L + i/2) (13) 

accounts for a regular behaviour near the origin for the modified 

J ' i( *Tl !° rr<< ! (14) 
\ nh{kr)) for r > 1 v ; 

In terms of a Fourier representation, using oscillator states as a basis, one then has the following well-defined 
expansion: 



7D 



\&W(r) — V" < r \ m> 

m=0 

oo 



This provides a proper definition for the coefficients cit^ and ci ' in (10), the "regular", resp. "irregular" 
asymptotic coefficients, the explicit form of which can be found in |2[ ffl. The equations for the asymptotic 
coefficients in Fourier space are then: 

oo 

Y < n | t - E | to > 4 +) = 0, 

m=0 

oo 

<n\f-E\m>c& = 0o <5«,o (16) 

m— 

One notices an identical form for all equations, except for the single one with n = 0. 

Substituting the solutions of ( |l6|) for the asymptotic coefficients in the set of equations ([[]), taking into 
account the substitutions @ and ([h]), then leads to the following linear system: 

oo oo 

^<n\H-E\m>c$ + tan(5)[#j <5„. + ^ < | V \ m > c^} = 

m=0 

oo 

< n | V | m > 4+» (17) 



m— m— 

oo 



m=0 



Equation ( pT[ ) now shows the influence of the potential energy matrix elements on the behaviour of the system 
in a clear way. 

Let us introduce the "Dynamical Regular and Irregular Coefficients" Vn and \Q as follows: 



oo 

= ^ < n | 1/ | to > cH (18) 



m=0 



5 



or, in an alternative integral representation in three-dimensional coordinate space: 

V( +) = I dr r 2 $ n (r) F(r) * (+) (fcr) 
= I dr r 2 $„(r) V(r) ^\kr) 



where we use the coordinate representation of $„(r), the n-quanta oscillator function < r \ n >: 



$„(r) = (-1)" N nLP L cM-\p 2 ) Ln + Hp 2 ) 



T 

(P=~ b ) 



(19) 



(20) 



2 (n!) 



1 



) 



]] T(n + L + 3/2) b 3 / 2 
Substitution of the Dynamical Coefficients (|l^) in ( |l7| ) leads to 

OO 

<n\H-E\m><$+ tan(<5) \/3 S n>0 + 



(21) 



m=0 



No approximations have been made so far to obtain this representation of the AM dynamical equations. 
We have taken into account the main (asymptotic) behaviour of the solutions in the equations, and used a 
regularization scheme to achieve this in a well-defined way. 

To reveal how the AM linear system can be solved in a numerically optimal, or at least acceptable, way, an 
analysis of the Dynamical Coeeficients Vn and Vn ' seems imperative. It is indeed clear that, if e.g. Vn + ^ and 
Vn \ from a given index n on, are sufficiently zero when compared to the other terms in (|2l]), the equations 
are reduced in a controlled and secure way. 

There is another important reason to investigate the behaviour of V^ + ' . It is indeed well-known from the 
integral equation formulation of quantum mechanics fill that: 



tam5r 



mk 



dr^^+\kr)V(r)^(k,r) 



(22) 



where ^(k,r) is the exact solution of the Schrodinger equation obtained with potential V. In an oscillator 
representation, (p||) reads as: 

. oc 

mk 



tan Sl = 



71=0 

r(+) 



(+), 



(23) 



from which we immediately recognize the importance of Vn in the convergence of the solutions of the contin- 
uous spectrum. 

The analysis of Vn + \ Vn ' must be done in terms of the basis index n, of the oscillator parameter b fixing 
the basis, and of the potential energy parameters. Results emerging from such an analysis can then be checked 
by studying the solutions ch^ and tan(<5) as a function of the same parameters. 



4.2 A Model Analysis using the Gaussian Potential 

To analyze the behaviour of the AM set of linear equations and its defining quantities in a more or less general 
way, we consider a simple gaussian potential of the form 

V(r) = V exp(-(-) 2 ) (24) 
a 

There are a number of reasons to justify such a choice: (1) an operator with gaussian functional form is easily 
handled in a harmonic oscillator basis, as matrix elements can be calculated using closed expressions, or simple 
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recurrence formulae; (2) (semi-) realistic potentials are often expressed explicitly as a finite sum of gaussians, 
each with specific amplitude and width parameters, and are of common use e.g. in nuclear physics calculations; 
(3) a very large class of potentials can be expressed in terms of a gaussian transform, such as e.g. a Yukawa or 
a Coulomb potential. 

In general the matrix elements of the potential energy operator, due to the r/a dependence of the latter, 
will depend on the ratio b/a. In the specific case of a gaussian potential one has an explicit dependence on 

7=A 2 (25) 
a 

This means that matrix elements of the gaussian potential energy operator in an oscillator basis are invariant 
under all scale transformations of a and b which preserve 7, and different physical situations will lead to the 
same matrix < n \ V \ m >. In particular, small values of 7 are realised by small values of b relative to o, 
or large values of a relative to 6; this situation corresponds to a "long range" potential. Large values of 7 
correspond likewise to a "short range" potential. 

The potential is certainly not the only parametrized quantity characterizing the set of AM equations. The 
oscillator radius of the oscillator basis is another and equally important parameter, because of the repercussion 
on the values of both the kinetic and potential matrix elements appearing in the equations. In this section we 
will therefore study the behaviour of the AM equations as a function of both the oscillator radius of the basis, 
and the width of the gaussian potential. Actually, as indicated by ( p5| ) only the ratio of these quantities is 
required, without any loss of generality for the current analysis. We therefore make a specific choice of potential 
parameters, namely Vq = —8 MeV and a = 1 fm, so that the potential has both a discrete and continuous 
spectrum. The 7 ratio will then be varied by adapting the oscillator radius b. 

A central theme in the analysis is the characterization for a rapidly converging, and stabilized, solution of 
the AM equations. This is most naturally done by searching for a maximal, limiting, value for the number 
of basis states involved in a specific solution, depending on the problem parameters used. As the asymptotic 
behaviour of the solution is governed essentially by the kinetic energy term, one should therefore take into 
account how both potential and kinetic matrix elements behave with respect to one another. In other words, 
it is not sufficient to know about the insignificance of potential matrix elements in absolute terms to decide 
whether one has reached the asymptotic region, but rather consider some relative insignificance with respect to 
the value of the kinetic matrix elements. 

When working in a coordinate representation, one can obtain a lot of information concerning the wave- 
function by analyzing the Hamiltonian only, and this both in terms of coordinate and energy ranges. The main 
reason for this is that the Hamiltonian is globally defined in the whole coordinate space. In a L 2 representation, in 
our case using an oscillator basis, this is not such a straightforward matter. The main parameters determining the 
rate of convergence, as well as the behaviour of the solution (the wave- function in an oscillator representation), 
are of a "global" nature, such as ^ m < n\V \ m > c m = < n \ V \ ^ >, < n \ V | > and so on. The 
matrix elements of the kinetic, respectively potential energy operators are "local" quantities (cfr. the value 
of these operators in a single point in coordinate space). The study of these elements however reveals the 
peculiarities of the "global" quantities, and will help to understand their behaviour. 

Concerning the kinetic energy term, two remarks are very important in the context of an oscillator basis. 
The kinetic energy contribution has a very simple representation in terms of the oscillator basis: (1) it is 
proportional to the inverse square of the oscillator radius &, and (2) the kinetic energy matrix has a tri-diagonal 
form, i.e. has non-zero matrix elements only along the major diagonal and the first subdiagonals. 

4.2.1 Analysis of the Hamiltonian matrix elements 

Fig. [I] displays the overall (qualitative) behaviour of a typical gaussian potential matrix in some specific 
oscillator basis. The main properties to be noted are the comportment of (1) the main diagonal which falls 
off monotonically for large n, and of (2) the non-diagonal matrix elements with non-zero values concentrated 
symmetrically around the main diagonal. This structure of the potential contribution is certainly of a promising 
nature, regarding the earlier remarks concerning the kinetic energy contribution. 

A more quantitative view of the diagonal behaviour of the potential matrix is displayed in Fig. ||, and this for 
a number of choices of oscillator bases. Three different views are presented in this figure: (a) the pure diagonal 
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Figure 1: Matrix elements of a gaussian potential in an oscillator basis with width parameter b — 0.75. The 
horizontal axes are labelled by the basis index, the vertical axis is in MeV. 

matrix elements; (b) the diagonal matrix elements multiplied by the square of the oscillator radius (eliminating 
the kinetic energy dependence on the choice of basis), and (c) the diagonal matrix elements normalized with 
respect to the first one; both figures || (a) and || (b) lead to Fig. || (c) through this normalization. From these 
figures it should be clear that one has to be careful when drawing conclusions for a proper choice of basis. Fig. 
U (a) might suggest that very large values of the oscillator radius b are optimal. Fig. || (b) on the contrary might 
suggest that very small values of b are optimal. The normalized Fig. || (c) finally suggests that, for the current 
potential parameters, the interval for b between 0.25 fm and 2.0 fm would be optimal. Fig. || at least reveals 
that there are important differences in behaviour of the potential energy contribution for different choices of 
the oscillator radius. Later on it will be shown in a more quantitative way that the qualitative conclusion for 
an optimal b interval as suggested by Fig. ^ (c) is correct. An interesting conclusion that is apparent from 
the figure is that a (coordinate space defined) short-range potential turns out to have a long-range character in 
an oscillator representation. As the oscillator basis is used as a Fourier basis to portray the solution, this is a 
well-known effect in terms of Fourier representation theory. 

A quantitative view of the non-diagonal behaviour of the potential matrix for the same bases displayed in 
Fig. H can be found in Fig. |^, and this by showing a typical row for a fixed column index (n — 50 in this 
case). The same 3 views (pure, multiplied by 6 2 , and normalized with respect to the diagonal matrix element 
with n = 50) are shown respectively in figures [| (a), (b) and (c). Again a strong dependence on the choice of 
the oscillator radius is remarked, but it is much more difficult to draw conclusions for an optimally converging 
basis from these figures. As this figure indicates that the potential energy contribution is concentrated around 
the main diagonal of the matrix, but with a relatively important distribution, it should already warn against 
carelessly applying the simple but straightforward solution scheme presented in section 2 of this paper! 

The previous figures thus only allow us to initiate some general discussion, but not to draw final and well- 
formed conclusions on the convergence problem, let alone on the determination of an optimal basis and the 
number of states involved in a stable solution. 
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Figure 2: Diagonal matrix elements of a gaussian potential in oscillator bases with different width parameter 
b: (a) exact values in MeV, (b) values weighted with b 2 in Mev fm 2 and (c) values relative to < | V \ >. 
The horizontal axis labels the basis index. 
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Figure 3: Matrix elements of a gaussian potential in oscillator bases with different width parameter b for row 
index n = 50: (a) exact values in MeV, (b) values weighted with b 2 in Mev fm 2 and (c) values relative to 
< 50 | V | 50 >. The horizontal axis labels the (column) basis index. 
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4.2.2 Analysis of the Dynamical Coefficients 



Up to now only local views of the potential contribution, i.e. matrix elements and their behaviour, have been 
considered to provide possible hints for a properly converging solution. The new representation of the AM 
equations presented earlier in this section introduced some new, global quantities which might be of interest in 
the convergence analysis. These are the Dynamical Coefficients Vn + ^ and Vn , and they combine all matrix 
elements of row n as indicated above. If, and when, these quantities would be (sufficiently) zero from a given 
value N on, this would certainly determine the maximal number of oscillator states to be considered for a 
converging solution. As Vn and V 7 [ depend on both the potential and the basis parameters, they will be 
analyzed as a function of the ratio of potential width to oscillator radius. As Vn + ^ and Vn~^ also depend 
explicitly on energy, this will be an additional parameter to consider. 

A qualitative view of Vn is shown in Fig. (a) as a function of both the index n and the ratio b/a; 



the quantity was normalized with respect to Vq , and the scale is a logarithmic one. From this figure one 
immediately notices that a well-defined value of the ratio of the potential to oscillator radius zeroes the quantities 
for already very small values of n. Fig. ^ (b) depicts the same view for Vn , carrying an analogous conclusion. 
Although both optimal values do not coincide completely, they are sufficiently close to define a narrow interval 
of optimal values for a swiftly converging solution. Indeed, if we consider a value N for which | /Vq \ < e 

and | V^/Vq ' | < e (with e.g. e < 10 -6 ), we can assume to have reached the asymptotic region. The value N 
then determines the number of basis functions needed to obtain a well-converged solution. Indeed, the solutions 
of ( |l7| ) are deviations with respect to the (known) asymptotic solutions, and become zero when both Vn and 
Vn are (sufficiently) zero. In other words, the value N from which on vi + ' and Vn , and thus the solutions 
Cn \ are zero determines what would be called in RGM terminology the matching point between the internal 
and the asymptotic region. It is important to remark that N is determined prior to solving the AM equations, 
and is obtained from the simple knowledge of the potential energy matrix elements. 

To corroborate the fact for an optimal choice of basis for a given gaussian potential, a quantitative view of 
and Vt ] as a function of the radius parameter b is shown in figures || (a) and (b) for a selected number of 
n values. From this figure one can determine for which values of the radius parameter one obtains a properly 
converged solution, by considering the intersection of both optimal intervals for Vn' and Vn . One notices 
that for a potential consisting of a single gaussian term, a very limited number of basis states is necessary for a 
proper convergence. For more intricate potential forms this will not necessarily be the case, although an optimal 
value for the radius parameter, be it associated with a relatively larger value of N, will still be available. 

To check the dependence of the so-called "optimal" b values (or regions) as a function of energy, Fig. || 

displays the values of both Vn + ^ and Vn~\ relative to respectively V^ and Vq~^ for varying b and energy for 
a fixed n. It is important to notice that, at least for the gaussian potential, the optimal value for the oscillator 
radius is independent of energy. 

Finally Fig. |?] shows the quantitative behaviour of Vn and Vn , again relative to respectively Vq and 

Vq \ as a function of n. From these figures one easily obtains an optimal number of basis functions, given the 
specific parameters of the problem and the precision considered. As a hands on example, Fig. |?] indicates that, 
for a precision of about 10 -8 , less than 12 basis functions are needed for all b values between 0.5 fm and 1 fm. 

The conclusions above indicate the importance of the Vn and Vn ' quantities, and a closer investigation 
imposes itself. A closed expression for can be obtained in a straightforward way as: 



(26) 



for which the full calculation is reproduced in appendix A. 

This expression features a single minimum in terms of 7, namely for 7 = 1/2 or b = a/y/2, with a value 

V M = y oiVni iI(^^+V2 exp( _l fc2) (2?) 



11 



9 



1.25 
1 . 5 



b (fm) 




0123456789 10 

b (fm) 



0.6 



0.5 







[ 


T 

1 


I I 1 1 1 


i i 


1 1 

n = 


5 




1 






n = 


10 




1 

1 




(b) 


n = 
n = 


15 - - ■ _ 

20 




1 
1 

1 
1 






n = 


25 




1 

1 
1 
1 

,1 
,1 
i| 
'1 












'1 
',1 
,1 
,i| 

[.M 




- 1 1 


i i 















0123456789 10 

b (fm) 



Figure 5: Vw (figure a, top) and V r i (figure b, bottom) as a function of b for various values of n and an 
energy of 1 MeV. 
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Figure 7: Logarithmic plot of 



K (+) /^o (+) | 



(figure b, bottom) as a function 



(figure a, top) and 

of n for several values of the oscillator radius b and an energy of 1 MeV. The vertical axis shows orders of 
magnitude. 
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which drops off to zero very fast in terms of n because of its following behaviour: 



1 

£+1/2 



(28) 



Both position and behaviour of V r [ + ^ around this minimum were already apparent from the foregoing figures. 
The asymptotic behaviour of V„ for both small and large values of 7 is also readily obtained. For small 
values of 7 (i.e. small values of the oscillator radius b compared to the potential width a) equation (|6|) yields 
the following asymptotic form, valid for large n: 



(+) 



VW « V exp{- 7 i?2 } d 

« Vo exp{- 7 i?2} yfm^j^kbRn) 



(29) 



This expression consists of two factors in coordinate representation, the gaussian potential V(r) and the (asymp- 
totic) wave function r), both evaluated in one and the same discrete point r = bR n (R n is the classical 
turning point, cfr. (||) and (0)). For 7 approaching zero the expression shows a slowly decreasing behaviour of 
Vn as a function of n. 

For large values of 7 one obtains the following asymptotic form valid for large n 



R-, 



17^} ^iM^Rn) 



(30) 



which for the limiting case of 7 now approaching infinity again displays a slowly decreasing behaviour as a 
function of n. 

Both very large and very small values of 7 will thus lead to rather badly converging solutions, which 
essentially means that a large number of basis functions will have to be considered when 7 reaches towards 
limiting values. 

This discussion on asymptotic behaviour confirms conclusions already apparent from the numerically calcu- 
lated figures of vi + \ 

Fig. H indicates that the asymptotic formulae are quite valid for a broad range of b values, and as such, can 
often be used to obtain reasonable values for N for properly converged solutions. Indeed, as the main factor 
which defines the decreasing behaviour of Vn for small, respectively large values of 7, is the potential term, 



which has the value exp{— jR n }, respectively exp{— ^}, as can be seen in equations (|2 
considers e.g. 10 -6 to be a reasonable measure for precision, one can write approximately 



and 



So, if one 



V, 



(+) 



(+) 
) 

(+) 



Vr 



( + ) 



exp{- 7 i?^} = l(T 6 (small 7) 



r>2 

exp{ ^} = 10- 6 (large 7) 

47 



(31) 



which immediately leads to the respective values N sw 14/(4 7) (for small 7) and N m M7) (for large 7); for a 
value of b/a = 0.2 one obtains N » 90, and for b/a = 3.0 N » 125. 

The asymptotic behaviour of can be used to an even better extent, by taking it into account when 
solving the AM system of equations. This will be pursued in full detail in a forthcoming section. 

Although more intricate to develop, an analogous analysis can be made for Vn , leading to corresponding 
conclusions. 



4.2.3 Analysis of the phase shifts 

In Fig. H we show the phase shifts obtained at an energy of E — 1 MeV for different values of 7 (actually in 
all further results a = 1.0 /to, so that b and 7 coincide), as a function of the number of basis states involved 
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Figure 8: 



compared to the asymptotic formulae discussed in the text, as a function of n in a 



logarithmic plot. The vertical axis shows orders of magnitude 



in the calculation. One notices an important gain in convergence speed (and thus precision) when using the 
reformulated version of the AM equations ( pl| ) , although small and large b (or 7) values still require an important 
number of basis states. 

On the same figure the normalized Cn^ solutions are shown for different values of b (7), for a calculation 
involving 100 basis states. These results are normalized with respect to the norm W 

V" I r (0) I 2 

W = ^ n =° (32) 

En= \Cn\l 

which is also include in the figures. One notices from this value that in the "optimal b (7)" region very few, 
very small Cn^ contribute to the solution. 

The results displayed on Fig. |^ confirm our suggestion made by analysing the behaviour of V^ + ' and V^~\ 
i.e. that for optimal values of b (7) convergence is achieved with less than 12 basis functions. They also confirm 
that for b — 0.2 one can use approximately 75, and for b = 3.0 approximately 100 basis states for reasonable 
convergence; the latter were indeed overestimated by the numbers obtained from equation (|3l|). 

To corroborate the fact that the reformulated AM equations (|2l]) provide faster converging solutions, we 
show in Fig. [To] both the c„ and coefficients obtained in a calculation at E = 1 MeV and b = 3.0, with a 
total number of 100 basis states. This picture confirms the fact that the true solution indeed deviates only by 
a small amount from the asymptotic one. 

To show the dependence of the results on energy, we show in Fig. [ll] the phase shift obtained with 5, 10 
and 15 basis states compared to the exact results, and this with the original ("simple") version (^) and the 
reformulated version (^l|) of the AM equations. For b values deviating reasonably from the optimal value, the 
latter form of the equations is seen to be highly superior to the original one. For very small and very large b 
values, the convergence is still problematic. 
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Figure 9: Left side figures: phase shift S (degrees) as a function of the number of basis states involved in 
the calculations; the energy was fixed at 1 MeV. Dotted lines refer to the solutions obtained in the original 
("simple") formulation, full lines refer to solutions obtained in the reformulated version. Right side figures: 
normalized c^ 1 values (see text) for a calculation with 100 basis states. The norm W is given as a percentage. 
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n 

Figure 10: c n (dotted line) and cl°^ (solid line) coefficients obtained in a calculation involving 100 basis states, 
with b — 3.0 and an energy of 1 MeV. 
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Figure 11: Phase shift <5 (degrees) as a function of energy, obtained with 5 (dotted line), 10 (dashed line) and 
15 (dash-dotted line) basis states. The solid line corresponds to the exact phase shift. Left side figures: Phase 
shifts calculated with the original version of the dynamical equations. Right side figures: Phase shifts calculated 
with the reformulated version of the dynamical equations. 
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5 Stable Solutions for the AM Equations 



In the previous section we obtained interesting results concerning the convergence behaviour of the solutions of 
the AM equations in terms of the parameters defining the basis and the potential energy. It was shown for a 
gaussian potential that the rate of convergence was reliably predictable. 

The parameter 7 or, for fixed potential width a, the oscillator radius b determines the rate of convergence 
of the expansion \P = J2 n c n \ n >■ It was indeed seen that an optimal value (or a restricted range of optimal 
values) for b leads to very fast convergence. However, values deviating more or less strongly from the optimal 
b lead to slowly converging and numerically unprccise results. 

As it is not always possible in a realistic calculation to choose the optimal b, for physical as well as for 
numerical reasons, it is important to develop strategies for stable results even in non-optimal 6-regions. This 
problem was already recognized in Q, were a coupled channels calculated for i He was performed with the AM, 
in which both cluster and collective configurations were taken into account: the b value was fixed by physical 
arguments (essentially optimizing the cluster channel results), and very non-optimal for the collective channels 
(in particular the Monopole Mode). 

The following subsections indicate how such problems can be solved by considering novel strategies to solve 
the AM equations, without modifying the original physical problem (e.g. by tampering with b). 

5.1 A General Analysis of the Asymptotics 

We attempt to evaluate the asymptotic behaviour (i.e. for n ^> 1) for the expansion coefficients c n =< n \ ^ > 
and the matrix elements < n \ V \ 'J >, using only very general information on the (unknown) wave function 
^. Based on this formulation we will then introduce strategies to overcome the slow-convergence situations. All 
assumptions made, and results obtained, will be checked against those produced with a gaussian potential. 

To derive a general view on the asymptotic behaviour of the Fourier coefficients c„, we use a Generator 
Coordinate (GC) representation for the wave function 4": 



using a scaled gaussian < r | (3 >= r L exp{— /3 2 r 2 } as a kernel of the integral transformation. It is well known 
that this kernel provides an alternative, continuous, basis to the traditional, discrete, oscillator basis. A tacit 
assumption in ( |33| ) is that the continuous spectrum wave functions are representable by such a representation, 
(in particular, that a "reasonable" form for g((3) exists, leading to an integrable result). 

5.1.1 Asymptotics of c„ 

The gaussian kernel of ([53|) is easily expanded in terms of the oscillator basis used throughout this paper. The 
basis functions | n > are explicited in a coordinate representation by (p0|). 

By substituting pb for r in the integral kernel r L exp{— (3 2 r 2 }, one can recognizes the generating function 
for the oscillator functions: 




(33) 




(34) 



where 



I3 2 b 2 



1 1 - £ 
2 1 + e 



or 



e = 



1 - 2f3 2 b 2 
1 + 2f3 2 b 2 



(35) 



The expansion coefficients < n | (3 > of exp{— /3 2 r 2 } in the oscillator basis are then 



<n|/3> = (l + e) L+ ^ 2 ^- b L 



n 



(1 + 2/3 2 & 2 )™+ L + 3 /2 N n 



(1 - 2(3 2 b 2 ) n 2 L + 3 / 2 foL 



(36) 
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from which one obtains the expansion coefficients c n =< n | ^ |> of ^ in the oscillator basis as 

< n | * >= / dj3g(f3) <n\[3> 



(1 - 23 2 b 2 ) n 2 L+3 / 2 

**iiTP^iinr <37) 

If, and when, <?(/?) is concentrated in a small vicinity of j3 w /?0j one can expect a highly convergent expansion 
for an oscillator length 6 w l/y2/3o. 

To study the asymptotic behaviour of the expansion coefficients (|37|), we consider two limiting regions for &, 
i.e. "small 6" and "large 6" . To this end, we rewrite ([37]) as follows 



where O = l/(V2b). 



fPo poo 

<n|^>=/ dflg{p) <n\p> +/ rf/3 3 (/3) < n \ (5 > (38) 



Small b values: for < < l/(y/2b) and large values of n, one can use the following approximate formula 

(1 ~ 2/32&2) " ~ exp{-i? 2 /3 2 & 2 } (39) 



(1 + 2/3 2 fe 2 )"+ i + 3 / 2 
For large values of n (n> 1) one also has 



iV» V 2I > + !) 

where again i?„ is the classical turning point (cfr. eqs (||) and (Q)). 

For very small values of b, one can consider only the first term in ([38|); using the approximations (|3^) and 
( p0[ ) one then obtains the following asymptotic form of ( p7| ) 

POO 

< n | >« v^/R^ / d/3<?(/3) (i?„6) L exp{-i? 2 /3V} (41) 
Jo 

which, considering (|33|), leads to the final approximation: 

c„ = < n | * > « vV^*^ i?„) (42) 

For small values of the oscillator length b, the expansion coefficients of the wave function \I/(r) on an oscillator 
basis are thus proportional to the wave function itself taken at the discrete argument values bR n . 

Such a relationship between the wave function in coordinate space and its oscillator representation was 
already obtained long ago [Q, Q for the regular and irregular asymptotic solutions. Later, this correspondence 
has been used as a heuristic principle for solving e.g. the Coulomb problem in an oscillator representation 

Large b values: for l/(y26) < /? < oo and large values of n, one can use the following approximate formula 

(l-2/3 2 6 2 )™ _ 1 (1 - l/2/3 2 6 2 )" 



(1 + 2/3 2 fo 2 )"+ i + 3 / 2 v ; (2/3 2 6 2 ) i+3 / 2 (l + l/2/3 2 6 2 )"+ L + 3 / 2 

« (-l)"7^^exp { -^ } (43) 

The expansion coefficients can then be approximated by the second term in ( j38| ) , which for very large values 
of b leads to the following asymptotic form of ( |37| ) 

< n | * > « (-i r V2X+V2 & x ^ ^ g ^ exp{-JSy> {u) 
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Using an integral transformation jl2| 11.4.29] the asymptotic form of the c n for large b becomes 

( Rr, 



Cn=<n\y> « (~l) n ^=^J^ drr 2 j L l^r)V(r) (45) 

= (-1)«1^(t) (46) 
V?r o 

where $ is the wave function in momentum representation. 

For large values of the oscillator length b, the expansion coefficients of the wave function ^f{r) on an oscillator 
basis are thus proportional to the wave function in momentum space taken at discrete arguments which are 
the momentum values R n /b. As oscillator functions in coordinate and momentum space only differ by a phase 
(— l) n and the argument p (which in coordinate space equals r/b and in momentum space equals kb), both wave 
functions ^(r) in coordinate representation and $(p) in momentum representation have essentially the same 
expansion coefficients. 

By integrating by parts, assuming that ^(r) is smooth and without singular points, and taking the 
leading asymptotic term for c n (i.e. for R n ^> 1), one obtains 

<n|*>« (_1)«-L*(0) (47) 

Rn 

The same result can be obtained by considering that the wave functions $>(p) of discrete and continuous 
spectrum states in momemtum space decrease at least as 1/p 2 - Substitution of this limiting factor in (^6|) again 
leads to © 

Summing up one notes that the asymptotics of c„ are related to the wave function behaviour in coordinate 
space: for small values of b the asymptotics are proportional to the wave function for large r, and for large 
values of b the asymptotics are proportional to the wave function in the vicinity of the origin. 

5.1.2 Asymptotics of < n | V | * > 

Asymptotic expressions for the matrix elements < n \ V \ $ > can be obtained in an analogous way as for 
c n =< n | "$! > by using the GC representation of the solution (]3^). We omit the details of the calculation, and 
only present the final results. 

Small b values: for small b values (or better, small 7 values), the matrix elements factorize as follows: 

< n I V I * > « V{bR n ) < n \ * > 

w V(bR n )V2^R n ~V(bR n ) (48) 

One notices that the behaviour of the matrix elements coincides with the one of the vE'(r), up to a trivial factor 
which is potential dependent. 

These results obtained for a general form of a short-range potential, are confirmed by the asymptotic form 
for Vn + ^ , evaluated earlier (^9|) for a gaussian potential. 

Large b values: for large values of b (or 7), the asymptotic form of the matrix elements < n \ V | \& > reduces 
to the integral 

< n I V I * >~ (-l) n ^^/Rn~ / drr 2 j L (^r)V(r)^(r) (49) 

Jo b 

or, in other words, to the convolution of the wave function with the potential in momentum representation 

<n \V \V >^ (-l) n ^^/R~ dkk 2 v(k n -k)<5>{k) (50) 
\/7r In 
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where v(k) and are the Fourier transforms of the potential V(r) and the wave function \I>(r) respectively. 

These results for a general potential can again be checked by considering the corresponding form for 
obtained earlier with a gaussian potential. Indeed, by calculating the integral 

V J drr 2 j L (-^r)exp {-jr 2 } j L (kr) (51) 

where the unknown solution \P(r) is substituted by \I>( + )(fc,r) = y/2/Trji,{kr) (the "free particle" solution), one 
obtains (^p|), which is valid for both large b and n. 

The matrix elements < n \ V \ ^ > can be evaluated in the asymptotic region by making the same 
assumptions concerning the behaviour of \&(r) as in the previous section to obtain the asymptotic form of 
c„ =< n | jfr >. Expanding \t(r) in a power series of r, keeping the first term ^(0) r L and integrating term by 
term in (|49|) leads to the following asymptotic form: 

< n | V | * >« (-l) n 4=V^(^r)*(0) (52) 

where the Fourier transform v(R n /b) of the potential determines how fast the matrix elements decrease in terms 
of n. 



5.2 New Strategies for Solving the AM Equations 

Based on the results obtained above, we can now suggest new strategies for solving the dynamical equations of 
the AM, in those cases where the oscillator length b, or more exactly the ratio 7 = (b/a) 2 , is relatively large or 
small compared to the optimal value. 

Small b values: we start from the Schrodinger equation in ([!]) 

J2<n\f-E\m> c m + < n | V | * >= (53) 

m 

^<i/|T , -S|^>c Al +<^|X>|*>=0 (54) 

where the index n is connected with the internal region and v with the asymptotic region, and N delimits the 
internal region. In the asymptotic region, we can use the asymptotic form (E3) for the matrix elements 



<v\V\*>= V(bR v )c v (55) 

which yields 

<n\f - E \ m> Cm+ < n | V | * >= (56) 

m 

Ys<v\T + V{bR v )- E\ fi>c ll =0 (57) 



In the original version of the AM the asymptotic coefficients c„ were given by (g), and obtained from a 
simple three term recurrence relation, due to the very selective coupling induced by the kinetic energy operator 
between oscillator states. Instead of the original form, and based on (|57|), we propose a modified three term 
recurrence form, including the asymptotic behaviour of the potential matrix elements in the diagonal term 

J2<"\T + V{bR v ) - E I fx > c M = (58) 
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or 



[<v\T\v> + V{bR„) - E]c u 

+ <v\t\v-l> c„_i + <v\t\u + \> c„+i = (59) 
In the asymptotic region c„ = c£ + tan<5c|> , hence the system (|59|) should be solved independently 



for regular (c£ ) and irregular (ci coefficients, with "boundary" conditions at starting remote points v = 
N a ,N a + l 



c« = ^2R v j L (kR v ) (60) 
c ( - } = y^2R^n L (kR,) (61) 

These modified asymptotic coefficients and ci ^ should be considered when calculating the correspondingly 
modified Vn and Vn in the dynamical equations (pl|). The modified equations should then be solved for the 



phase shift S and the coefficients of the internal region . 

Large b values: for large values of b, we consider the dynamical equations in the form (0), containing the 
vS +) and terms 



Y,<n\H-E\m>cW+Vt ) tan(<5) = - V« 



(62) 



for the internal region [n < AT), and 



2<i/|f'-.E|f«>c M +<i/|tf'|¥>=0 



(63) 



for v > N. 



Substituting effl + cj^ + tan(5)c)j 1 for in (63), an assuming that cjjf' and J satisfy the three-term 



„(-) 



( + ) 



„(-) 



recurrence relation (|T^), one obtains 



^2<v\f-E\n> < 1/ 1 V" | * >= 



(64) 



Taking the asymptotic form ( p2[) into account, one obtains the following equation for c 



(a) 



^<z,if-£;| A1 > c (' i ) + (-ir v ^: 



^=0 



(65) 



where c£, was substituted by 



(66) 

introducing the " asymptotic" coefficients and the " residual" coefficients ci . All constants in the asymptotic 
form (52) were omitted and replaced by the global factor W to be determined. Equation (^35|) should now be 
solved subject to the boundary condition 

ci a} =0 (67) 

for v = N a , N a + 1 . Having obtained , the equations for the internal region should be modified to take the 
substitution (pq) into account 



<n\H - E\m> c$ 



m=0 



y(<0 _ Ec (a) 



W + tan((5) = - V^ +) 



(68) 
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were vS a ^ now stands for the sum 



K (a) = E < n I # I m > c m ( 69 ) 



m— 



The meaning of the asymptotic coefficients c£ can be most easily understood by considering the major 
part of the asymptotic behaviour already in an intermediate n region. Indeed, for large b values, there is a 
very slow convergence of the results. This implies that the internal region extends towards very large n, with 
corresponding cl°^ solutions probably very close (i.e. to first order) to the el . 

Numerical application: To demonstrate how the new ("small b" and "large 6") strategies accelerate the 
convergence, we again consider a gaussian potential in the two limiting cases where (i) the oscillator radius is a 
factor of 5 less than the potential range, and (ii) the oscillator radius is a factor of 5 larger than the radius of 
the potential. 



Fig. 12 compares the phase shifts, obtained at an energy of 1 MeV, obtained in the original, reformulated 
and asymptotic approaches. 

In Fig. |l^ we display the exact phase shifts and those obtained with the new strategies using 5, respectively 
10, basis states. One notices that for the "small b" case, 5 basis functions yield almost exact results (within 
a precision less than 0.01%); by even considering only one basis function, the phase shift is obtained within a 
precision of 1%! In the original form the latter precision could only be reached by using more than 50 basis 
functions. For the "large b strategy" (b/a = 5), the convergence speed is remarkably increased, although less 
spectacular than in the "small b" case. The number of states needed for an identical result in the original 
formulation (|2l|) is about 5 times higher. 

The two strategies suggested are thus seen to significantly improve the convergence of the results and to 
reduce the computational efforts to obtain a desired precision. 

By modifying the form of the asymptotic AM equations through inclusion of dynamical features as proposed 
and realized above for both small and large values of b, one is naturally led to the introduction of an "intermediate 
region" . This region distinguishes itself (i) from the "internal region" where the solutions are governed by the 
potential, and (ii) from the "asymptotic region" where the kinetic energy dominates the equations. In the 
intermediate region, the solutions are easily (i.e. in a numerically simple way) obtained. There are then actually 
two variational parameters to consider, N marking the border between the internal and intermediate region, 
and N a marking the transition from the intermediate to the true asymptotic region. The larger or smaller is 
7, the larger is the size (i.e. the difference N a — N) of the intermediate region. The strategies proposed above 
for small and large values of 7 have shown that, by considering a sufficiently large intermediate region (e.g. 
N a = 500) , significantly reduces the range of the internal region, and dramatically decreases the computational 
effort to obtain converged results. 

The intermediate region is charactarised by a noticeable but non-dominant presence of dynamical (i.e. 
potential) effects. These effects can be well approximated by asymptotic forms of both the expansion coefficients 
c n and potential matrix elements < n \ V \ ^> >. For small 7 (long-range correlations) , this leads to a redefinition 
of the equations of cl + ^ and cl ' by incorporation a potential term in the three term recurrence relation governing 
the true asymptotic behaviour. For large 7 (short-range correlations) one approximates part of the solution 
by its asymptotic behaviour ■ The latter are solutions of an inhomogeneous set of linear equations containing 
the Fourier transform of the potential, and as such accumulate the potential effects in the intermediate region. 

In the region of 7 near to the optimal value, where fast converging solutions are obtained, there is no need 
for an intermediate region, and thus for an approximation of the dynamical equations; the values of N and N a 
coincide in this case. 



6 Conclusions 

By introducing the "Dynamical Coefficients" Vn^ 1 and Vn ' we were able to suggest a new form of the dynamical 
equations of the Algebraic Version of the Resonating Group Method (AM). The dynamical coefficients allow one 
to determine oscillator basis parameters (oscillator length and number of basis functions) for obtaining otimally 



2G 




Figure 12: Phase shifts at an energy of 1 MeV as a function of the number of basis states involved in the 
calculation. Dotted lines represent solutions of the original ( "simple" ) formulation, solid lines represent solutions 
of the reformulated version and dashed lines represent solutions of the "new strategies" . 
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Figure 13: Phase shifts in the energy range from to 10 MeV for a ratio b/a — 0.2, obtained in the "small b 
strategy" (top), and b/a = 5, obtained in the "large b strategy" (bottom). Full lines show the exact solution, 
dotted (resp. dashed) solutions are obtained with 5 (resp. 10) basis states. 
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converging and stable solutions of these equations. It was shown that for an optimal parameter set 5 to 10 
functions were sufficient to obtain results with a precision of more than 99.9%. 

The new form of equations have also been shown to be a proper starting point for obtaining new solution 
strategies in parameter regions were slow convergence would occur. Indeed, in parameter regions were the 
oscillator length is much larger, respectively much smaller, than the optimal one asymptotic approximations 
could be formulated by analysing the AM equations. The latter can be easily implemented numerically, and 
have shown to dramatically improve the convergence of the solutions. 

It was shown also that, if the oscillator radius is much smaller than the width of the potential, the expansion 
coefficients c n of the wave function on the oscillator basis coincide with the wave function in coordinate space 
up to a simple factor. In the other limiting case, when the oscillator radius is much larger than the width of 
the potential, the expansion coefficients c n were found to be proportional to the wave function in momentum 
space. 
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A appendix 

In this appendix the explicit form, and governing recurrence relations, of the dynamical coefficients Vn +S> and 
Vn ^ are calculated for a gaussian potential. 
The following integrals should be evaluated: 




(70) 



(71) 



where <Ia + ) 



f j L (kr) and #(~) 



2 



is a solution of the inhomogeneous differential equation: 



(f =/3 $ (r) 
(f3o is defined in eq. ( |l3|) ), which can be represented in the integral form 



(72) 





(73) 



Integral (^TJ) 




l-4 7 2 



) 



(74) 
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If 7 = s, then 

V; (+) = ^L^(^) 2n+i+1/2 exp(-ifc 2 ) 
The second integral (|7l]), reduced to the form 



(75) 



can be expressed through the c„ ' coefhcients 



(1 + 27)™+( i+3 )/ 2 (l + 7 )™+ i / 2 

_n»-I + :V2) N nL 
ml". 

where 



" '-'/ffi (78) 



New variables fc' and k' 



*'-Vr^*- i '-Vr^ i - (79) 

and the following relation for Laguerre polynomials (see fTij , vol. 2, formula 10.12(40)) 

* — ' to!1 n — m + a + 1) 

m— 7 

were introduced in ( |73| ) to obtain this result. 

To obtain the recurrence relations for Vn + ^ and V„ \ we start from the equations for the functions 
and ^ (_) 

(T = 0, (81) 
(T -£)*(-> = A)*o(r) (82) 

Multiplying both sides of these equations by $> n (r)V(r) and integrating, one obtains 

< n\VT\¥ +) > -EV,i +) = 0, (83) 

< n\VT\^ > -EV^ = (3 <n\V\0> (84) 

On notices the similarity of the left hand side of both equations, leading to identical recurrence relations for 
both type of coefficients. The matrix element < n \ V | > for the gaussian potential is 



< n 



i v i o >= (-mm - «)-^ XAv^' < 85 > 



where 

z = ( 7 + l)- 1 (86) 
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Using the commutation relations of the operators T and V, or the explicit form of V„ , one finally obtains 

(1 + 2 1 ) 2 T rhn+1 v£\ + [(1 - 47 2 )T„,„ - E]V^ + (1 - 2 1 ) 2 T n ^ 1 V r i +\ 

= (87) 

(1 + 2 1 ) 2 T n , n+1 V^\ + [(1 - 4 7 2 )T„,„ - E]V.t ] + (1 - 2 1 ) 2 T n , l . 1 vt\ 

= Pa < n\V\0 > (88) 

If in ( pSj ) the potential is switched off (Vo = 1 and z — 1 (or 7 = 0)), the recurrence relations for V 7 { + ^ and 
Vn revert to these of c„ + ^ and c„~\ as is to be expected. 
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